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We calculate the yield of high energy neutrinos produced in astrophysical sources for arbitrary 
interaction depths to and magnetic field strengths B. We take into account energy loss processes 
like synchrotron radiation and diffusion of charged particles in turbulent magnetic fields as well as 
the scattering of secondaries on background photons and the direct production of charm neutrinos. 
Meson-photon interactions are simulated with an extended version of the SOPHIA model. Diffusion 
leads to an increased path-length before protons leave the source of size R s and therefore magnetized 
sources lose their transparency below the energy E ~ I0 18 eV (R s /pc) (B/mGjr 1 '", with a = 1/3 
and I for Kolmogorov and Bohm diffusion, respectively. Moreover, the neutrino flux is suppressed 
above the energy where synchrotron energy losses become important for charged particles. As a 
consequence, the energy spectrum and the flavor composition of neutrinos are strongly modified 
both at low and high energies even for sources with to < 1 . 
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I. INTRODUCTION 

High energy neutrinos from astrophysical sources are 
the decay products of secondary mesons produced by 
scattered high energy protons on background protons or 
photons. A classic example are the so-called cosmogenic 
or GZK neutrinos produced in scatterings of extragalac- 
tic ultra-high energy cosmic rays on cosmic microwave 
photons during propagation Since the sources of 

ultra-high energy cosmic rays unavoidably contain also 
matter and photons, all cosmic ray sources like active 
galactic nuclei (AGN) or gamma-ray bursts should pro- 
duce to some extent also high energy neutrinos 0, H, 0, H[ . 
However, the elusive nature of their interactions has 
prevented so far the detection of neutrinos from other 
sources than the Sun and SN1987A, despite of intensive 
experimental efforts. 

Theoretical limits or predictions for neutrino fluxes are 
therefore important as a guideline for experimentalists. 
Two different kinds of bounds on high energy neutrino 
fluxes exist: The cascade or EGRET limit uses bounds on 
the diffuse MeV-GeV photon background to limit the en- 
ergy transferred to electromagnetically interacting parti- 
cles that are produced unavoidably together with neutri- 
nos @. The cosmic ray upper bounds of, e.g., Rcfs. 0, [I[ 
use the observed ultra-high energy cosmic ray flux to re- 
strict possible neutrino fluxes. One of the several under- 
lying assumptions of the latter limit is that all neutrino 
sources are transparent to hadronic interactions and thus 
at least neutrons can escape from the source region with- 
out interactions. 

In a previous work, wc calculated the flux of high en- 
ergy neutrinos produced as secondaries in astrophysical 



sources with arbitrary interaction depth but negligible 
magnetic fields Q. The present work extends Ref. 
in two respects: First, multiple interactions that become 
relevant in a thick source require the modeling not only 
of photo- nucleon (N + 7 — ► X), but also of photo-meson 
(7r + 7 — > X, K + 7 — * X) interactions. The latter 
are treated now within a self-consistent extension of the 
SOPHIA model Second, and more importantly, we 
have added the effects of possible magnetic fields around 
the source region. This includes diffusion of charged par- 
ticles in turbulent magnetic fields and energy loss pro- 
cesses like synchrotron radiation. As a result, the shape 
of the neutrino spectra deviates strongly from the injec- 
tion spectrum of protons even for transparent sources, 
To < 1 , both at the low and the high energy end of the 
produced neutrino flux. Moreover, the relative impor- 
tance of the various channels contributing to the neu- 
trino yields changes strongly, leading thereby to large 
variations in the neutrino flavor composition as function 
of the energy. Specific examples are the cases where the 
neutrino spectrum is dominated by kaon decays [ll|, [l2j 
or influenced by muon damping [13l | . 

This work is structured as follows: In Sec. II, we dis- 
cuss how we simulate the interactions and the propaga- 
tion of charged particles. In particular, we describe our 
extension of the SOPHIA model, the continuous energy 
loss processes considered and how we account for diffu- 
sion. In Sec. Ill, we discuss the influence of magnetic 
fields on the escaping proton flux. Armed with our un- 
derstanding of the proton fluxes, we continue in Sec. IV 
to examine the resulting neutrino fluxes and in Sec. V 
their flavor composition. Finally, we summarize our re- 
sults in Sec. VI. 
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The application of our simulation to concrete astro- 
physical models and a comparison with earlier studies of 
neutrino fluxes from AGN cores Hand jets [H, clusters 
of galaxy [ijj], or hidden sources |15j will be performed 
in paper III of this series [l6| . 



II. SIMULATION OF PARTICLE 
INTERACTIONS AND PROPAGATION 



A. Simulation of photon-hadron interactions 

1. Photo-nucleon interaction in the SOPHIA model 
The SOPHIA model 

was designed to treat gamma- 
nucleon interactions in the broad energy range from the 
production threshold of mesons up to TeV center of mass 
(cm.) energies i/s. The underlying physics changes 
substantially over this energy range. At low energies, 
typically below y/s ~ 1 -j- 2 GeV, the interaction pro- 
cess is dominated by the excitation of resonant states of 
the incident nucleon and by the so-called direct meson 
production process. On the other hand, at higher ener- 
gies the process is dominated by the ^-channel exchange 
of composite states, Reggeons and Pomerons, leading to 
multiple production of secondary hadrons. 

The cross section for s-channel resonance production 
is described by the Breit-Wigner formula and results in 
a strong energy dependence with pronounced resonance 
peaks. All the basic parameters like resonance mass, 
width, branching ratios as well as angular distributions 
for resonance decays are taken from experimental data; 
when no information on the angular distributions was 
available, an isotropic decay was assumed. 

Direct meson production corresponds to a i-channel ex- 
change of a virtual meson state, which implies a strong 
coupling with the nucleon and an electromagnetic (e/m) 
one with the photon. As a consequence, neutral meson 
exchanges are strongly suppressed and one deals with 
charge exchange processes; dominant contributions in 
case of 7p collisions arc 7p — > 7r + n, jp — > 7r~ A ++ , the jn 
case follows from isospin symmetry. The contribution of 
heavier meson and nucleon exchanges are strongly sup- 
pressed as 1/mf, m t being the mass of the exchanged 
state, and can be neglected. Still, with increasing s ad- 
ditional contributions come from other ^-channel reso- 
nances lying on the Regge trajectory with pion quantum 
numbers, so that the corresponding contribution to the 
cross section behaves as s 2 ( q t(°)- 1 ) i where 0^(0) — is 
the intercept of the pion Regge trajectory. Therefore, in 
the SOPHIA model the contribution of the direct process 
to the 7^V cross section is described by a smooth func- 
tion, which vanishes near the hadron production thresh- 
old and behaves like s~ 2 in the high energy limit. The 
corresponding normalization was tuned to reproduce ex- 
perimental data on binary production processes. The 
angular distribution for the direct process is strongly for- 



ward peaked and was parameterized as 



i7rr dir 



dt 



(1) 



where t is the momentum transfer squared in the pro- 
cess. The value of the slope bdh- was taken from a fit to 
experimental data [l7j . 

At higher energies the interactions are dominated by 
Reggeon and Pomeron exchanges between the nucleon 
and the hadronic state of the photon, which result in mul- 
tiple production of secondary hadrons. Both Reggeons 
and Pomerons are composite exchanged states which 
emerge in the so-called 1/N C , 1/Nf QCD expansions [HI, 
where Reggeons correspond to a set of diagrams char- 
acterized by a planar topology, whereas Pomerons are 
described by non-planar (cylinder) diagrams. Qualita- 
tively, one can see a typical Rcggcon-dominatcd inter- 
action as follows. The incoming photon converts into 
a vector meson state, consisting of a "dressed" color- 
connected quark-antiquark pair. Next, the antiquark an- 
nihilates with a quark from the nucleon, such that both 
the spectator quark from the vector meson state and the 
spectator diquark from the nucleon are connected by a 
color field. While these spectator partons move apart 
with their initial velocities, the color field is stretched 
into a tube (color string) and finally disrupts, resulting 
thereby in the production of secondary hadrons, which is 
described by the string fragmentation model [19[ . On the 
other hand, Pomeron exchange does not lead to a trans- 
fer of quantum numbers between the interacting states. 
Instead, it induces a color exchange between the valence 
constituents of the vector meson and the nucleon, such 
that color strings are formed correspondingly between 
the quark from the meson and the diquark from the 
nucleon, and between the meson antiquark and the nu- 
cleon quark. Thus, Pomeron exchange results in the pro- 
duction and hadronization of two color neutral strings, 
compared to just one for the Reggeon exchange. Apart 
from the string fragmentation procedure, one has to spec- 
ify here the energy-momentum partition between valence 
constituents of the interacting hadronic states. They are 
defined by the intercepts of secondary Reggeon trajec- 
tories [20j and the distribution of light cone momentum 
fraction of valence quarks (x) and diquarks (1 — x) in the 
nucleon is approximated as 



1 



fq/N(x) ~ -= 



(1-x) 



1.5 



(2) 



while for the valence quarks (x) and antiquarks (1 — x) 
in mesons 



fq/-ir(p,u) ( x ) 



1 



\A(1 ™ X) 



(3) 



was used. 

Multiple Pomeron exchanges can be neglected still 
with reasonable accuracy in the SOPHIA energy range. 
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parameterized as [2c 



(a) 



(b) 



(c) 



FIG. 1: Direct tv + production in "/p (a) and jtt + (b), (c) 
interactions. 



However, the model takes into account diffraction pro- 
duction mechanism, which is treated phcnomcnologi- 
cally as a constant share of the Pomeron and Reggcon 
exchange contributions. The angular distribution for 
diffraction interactions is strongly forward peaked and 
is parameterized similarly to |T]), using an energy depen- 
dent slope frdiffr from pTj . 



2. Generalization for photo-meson interactions 

We have performed a direct generalization of the 
SOPHIA approach for the case of gamma-meson (tt , 
K^, Kq, Ko) interactions. As above in the 7-/V case, the 
dominant mechanisms are resonance production and bi- 
nary processes at low energies and Reggeon (Pomeron) 
exchanges at high ones. We have included the known res- 
onances, with well measured c/m branching width, from 
the compilation in the Review of Particle Properties [2l[ , 
the complete list is given in Table [IJ All the resonances 
were assumed to decay isotropically. Concerning the di- 
rect meson production, the only difference compared to 
the 7-ZV case is that the virtual meson-nucleon coupling 
is replaced by the meson- meson one, as shown in Fig. 1. 

In the energy range of interest, experimental data on 
binary processes in hadron-hadron and gamma-hadron 
collisions are in a very good agreement with the additive 
quark model (AQM) predictions [22|. Thus, we define the 
mcson-hadron coupling as the sum of the contributions 
corresponding to the coupling to "dressed" constituent 
quarks. Correspondingly, for the cross section for the 



binary process 777 



tt u we use the SOPHIA pa- 



rameterization for the reaction 775 — + ir + n (with prop- 
erly modified hadron production threshold), based on the 
fact that the number of valence (anti-) quarks in pion is 
the same as the number of valence u-quarks in proton. 
The same parameterization is used for r yK + — > ir + Ko 
(similarly for other kaons), now being rescaled by a fac- 
tor 1/2, corresponding to having just one light valence 
(anti-) quark in kaons. Here we neglect the f-channel ex- 
change of kaon states, which is mass squared suppressed 
compared to pion exchange, as discussed above. 

A remarkable property of Reggcon exchanges is the 
factorized form of the corresponding amplitude, which is 



q rC (*) 



1 + C exp(-i7r a R c (t)) 
sin(-i7ra H c (t)) 



(5) 



where a R c (t) is the Regge trajectory for the Reggeon K^, 
<7* (t) the Reggeon-hadron a vertex, and £ = ±1 is the 
Reggeon signature; so — 1 GeV 2 is the hadronic mass 
scale. This expression is also valid for Pomeron exchange, 
the latter being a particular Reggeon case, correspond- 
ing to the intercept ap(0) > 1, positive signature, and 
zero transfer of quantum numbers. This allowed to con- 
struct very successful parameterizations for high energy 
behavior of hadron-hadron and gamma-proton cross sec- 
tions [2J], which include only three terms: exchanges of 
Reggeons with positive (K + ) and negative (R - ) signa- 
ture, and of the Pomeron (P) 
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where the constants Cf b are related to the Reggeon 
and Pomeron couplings to hadrons a and b, C^ b ~ 
9a (0) 9b (0)' an d the energy exponents of the three terms 
in © are expressed via the intercepts ax(0) of the cor- 
responding Regge trajectories and do not depend on the 
types of interacting hadrons. Thus, we can obtain the 
corresponding parameters for 77r and jK interactions as 



a 



c. 
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(7) 

(8) 



for X = 



On the other hand, we have C* 



' 7 fc - 0. as 

the Reggeons with negative signature do not contribute 
to "fh scattering. The parameters C* p , , ax(0) have 
been taken from the same cross section fits [2l| which 
have been already employed in the SOPHIA model. 

The calculated total 77T and jK cross sections, as well 
as partial contributions of various interaction channels, 
are plotted in Figs. [2HH 

In contrast to the 7 A case [lfj, production of reso- 
nances does not play a very prominent role in 77r and 
collisions, because of the comparatively high energy 
thresholds of these processes. On the other hand, of 
considerable importance is the direct pion production, 
which starts already at the hadron production threshold. 
It is worth stressing again that both the direct contri- 
bution and the one of Reggeon and Pomeron exchanges, 
dominant at high energies, are well-defined by the proper 
rescaling from the 7A case; the former using the AQM 
picture, the latter within the Regge theory framework. 

For particle production, we have used as SOPHIA the 
Lund JETSET 7.4 string fragmentation procedure [l9l ]. 
as well as the distribution ([3]) for the energy-momentum 
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FIG. 2: The jit total cross section (solid line) and the con- 
tributions of resonance production (dashed line), of the direct 
process (dotted line), and of Reggeon and Pomeron exchanges 
(dot-dashed line). 



FIG. 4: The 7A0 total cross section and the contributions of 
various interaction channels. The abbreviations for the curves 
are the same as in Fig. [2] 
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FIG. 3: The jK ± total cross section and the contributions of 
various interaction channels. The abbreviations for the curves 
are the same as in Fig. [2] 

partition between valence constituents. Also we have 
kept the same proportion between diffractive and non- 
diffractivc final states for "fir and interactions as the 
one for case in SOPHIA 1 . In summary, the described 
generalization of the SOPHIA interaction mechanism for 
77r and jK collisions is well-defined: All the relevant 



TABLE I: The meson resonance processes considered and 
their characteristics: mass M, width V, e/m decay width 6 7 , 
and hadronic decay modes. 



resonance process 


M/GeV 


F/GeV 


10 3 6 7 


decay modes 


77T* - 


+ P±(770) 


0.768 


0.150 


0.45 


7T7T (100%) 


77T* - 


+ 6f (1235) 


1.230 


0.142 


1.6 


LOTT (100%) 


77T* - 


+ aj(1320) 


1.318 


0.107 


2.68 


rpr (14.5%) 












pTT (70.1%) 












LOTV (10.6%) 












KK (4.8% ) 


7Jf ± - 


-> A'* (±) (892) 


0.892 


0.0508 


0.99 


Kty (100%) 


7^ - 


^A 2 * (±) (1430) 


1.426 


0.0985 


2.4 


K-r (50%) 












A'*(892)tt 












(38%) 












Kp (9%) 












Ku (3%) 


7 A-°- 


■* AT* (0) (892) 


0.896 


0.0507 


2.3 


Kn (100%) 



parameters are fixed either by experimental data or by 
theoretical arguments. Finally, we note that in compari- 
son with the simplified modeling of pion and kaon inter- 
actions in Ref. [9(, the interaction lengths of pions and 
kaons are reduced outside the resonance region by ~ 35%. 

B. Interactions in the continuous energy loss 
approximation 



One may expect some difference in the diffraction excitation 
probability for pions and kaons compared to protons. The data 
on pp, up, and Kp interactions indicate that this difference is at 
the 10% level and therefore inessential for our applications. 



The following interactions are treated in the continu- 
ous energy loss approximation. 

1. Synchrotron radiation: Processes in external mag- 
netic fields are characterized by the parameters \ = 



5 



eF^ p Vp»/m 3 = (p±/m)(B/B a ) and B cx = m 2 /e. 
The intensity of synchrotron radiation can be ap- 
proximated for B <C B cr by (2f| 



W 

mi 



1 



55^ 



24 

2 8 / 3 r(2/3)r 4/3 



64 2 
3 5 



£ » 1 , (9) 



where £ = 2^/3 and W c ] = (2/3)am 2 x 2 is the clas- 
sical intensity of synchrotron radiation. We glue 
the two approximations at £ = 0.8 together. The 
energy loss per time is (3 = dE/dt = — W. We de- 
fine as typical length-scale l syn of synchrotron losses 
(in the classical limit £ <C 1) 



1 dE 
E ~dt 



3 1 



2am 7 



-B CI 



(10) 



2. Below the pion-production threshold, e + e~ pair 
production on photons becomes for low magnetic 
field strengths the main energy loss mechanism of 
charged hadrons. We use the energy losses as cal- 
culated in Ref. 

3. For muons, we include energy losses due to inverse 
Compton scattering from Ref. [271 ] . 

In the following, we will consider static sources and ne- 
glect the possible effects of adiabatic cooling as well as 
the energy losses due to curvature radiation. The latter 
is formally equivalent to synchrotron radiation after the 
replacement of the Larmor radius Rl = E/(eB) with 
the curvature radius Rq and needs therefore no special 
consideration. 



to specify thus the energy 
D(E). Following Ref. [2£ 



dependent diffusion coefficient 
, we use 



D{E) = D 



Rl 

17 



Rl 

17 



for E <E L , (13) 



where l c denotes the coherence length of the magnetic 
field and a = 1/3 and 1 correspond to the case of Kol- 
mogorov and Bohm diffusion, respectively [29| . 

The diffusion time idiff is defined as the typical time a 
charged particle needs to diffuse the distance R s assum- 
ing negligible energy losses, 



tdm - 6D 



(14) 



Thus the effective size R e g of a magnetized source be- 
comes i? ff = tdiff = Rs I (6-D) ■ We choose the normal- 
ization constant Dq in Eq. (TIH|) such that the effective 
source size equals the true one, i? c ff = R s , for El, 



Rs 



Dn = ^ 



lc. 



1 



6 \RJ 1 + {1 C /R s y 



(15) 



Without diffusion, one defines the interaction depth as 
the ratio tq = R s /k n t of the source size R s to the in- 
teraction length Z; nt of nuclcons with photons. For the 
illustration of our numerical results, wc determine in the 
following T via l- ult = l/(n 7 er) with a = 0.2 mb above 
the threshold -Eth as reference cross section. Therefore 
in the case of a background of thermal photons with 
temperature T, the interaction length becomes Zi n t ~ 
2.5 x 10 26 (T/K)~ 3 cm. In the case of diffusion, it is con- 
venient to introduce additionally an effective interaction 
depth r of f = Res/knt or 



C. Diffusion in magnetic fields 

The magnetic fields present in the surrounding of the 
acceleration region induce not only synchrotron radia- 
tion as an important energy loss process discussed in 
the previous subsection, but lead also to deviations from 
straight-line propagation for charged particles. The tran- 
sition between the ballistic and the diffusion regime hap- 
pens approximately at the energy El when the Larmor 
radius, 



R L = 1.08 x 10 _3 pc 



E G 



10 18 eV B ' 
equals the source size R s , i.e. at the energy 

E L = 10 18 cV (i? s /pc) (B/mG) . 



(11) 



(12) 



In a turbulent magnetic field, the path of a charged 
particle can be modeled either as a random-walk from a 
microscopic or as a diffusion process from a macroscopic 
point of view. We shall use the latter picture and need 



Toff = 



To 
T 



for E>E L 



(E L \ a i+{Rs/lcY- a fnr F , Fr 



(16) 



Thus a source that is transparent at high energies, when 
particles move in the ballistic regime, becomes thick in 
the diffusion regime below the energy ~ Elt^"". 



D. Summary of interactions and propagation 

We idealize a neutrino source as an acceleration region 
surrounded by a sphere of radius R s containing photons 
and turbulent magnetic fields. The photon energy distri- 
bution may follow either a (broken) power-law or a ther- 
mal Planck spectrum, although here we present results 
in most cases only for a Planck distribution. The prob- 
ability J\f that a particle diffuses outwards the distance 
Ar without scattering or decay is given by 



TV = exp 



r+Ar 



(17) 
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where l^ cc and Zi nt are its decay and interaction length, 
respectively 2 . The path length AZ of the trajectory and 
the distance Ar diffused outwards are connected by 



(18) 



Thus diffusion enhances interactions and energy losses 
by the factor u(E) oc 1/D(E) for E < E L . The energy 
of the particle along the path I is obtained by integrat- 
ing the energy losses (3(E). Note that in the calculation 
of diffusion coefficients one assumes generally negligible 
energy losses. Thus, Eq. (TT8"]) neglects the decrease of 
the diffusion coefficient along the particle trajectory and 
hence underestimates the interaction depth and the true 
number of interactions. 

In our Monte Carlo simulation, we track explicitly all 
secondaries (N, 7r ± , K^, if£ s ) for which the interac- 
tion rate is non-negligible compared to their decay rate, 
and take into account the possibility to produce prompt 
charm neutrinos; for details see @. 



III. PROTON FLUX 

For the illustration of our results we have chosen in 
the following always for the initial proton flux a power- 
law dependence, dN p /dE — KE~ ag with a g = 2, and a 
rather high value of the maximal energy E max = 10 24 eV. 
In order to facilitate the comparison with our previous 
work Ref. @, we consider only thermal distributions of 
background photons. 

Sources with negligible magnetic fields were thoroughly 
discussed in Ref. [9| . While in this case the initial and fi- 
nal proton fluxes deviate strongly for thick sources, they 
basically coincide for transparent ones. By contrast, the 
final proton spectrum is strongly distorted even for a 
transparent source, if the magnetic field in the source 
is sufficiently large: At high energies, synchrotron radia- 
tion losses lead to a change in the spectral index of the 
final proton flux, while at low enough energies protons 
start to diffuse in the turbulent magnetic field. Thereby 
the number of interactions increases, the proton flux is 
suppressed and at the same time the production of neu- 
trons is enhanced. As we will show now, the cosmic ray 
flux from a magnetized source has bumps and breaks, 
and may even suddenly decrease by orders of magnitude. 



A. Low-energy protons and diffusion 

We start our analysis by studying the consequences of 
diffusion on the final proton flux. In Fig. [5j we show the 
effective size R e g of a thin source with T = 10 4 K, R s = 
2.5 x 10 13 cm (r = 0.1) and a magnetic field strength 




The corresponding equation in Ref. [j| contains a typo. 



E(eV) 



FIG. 5: Effective size (thick solid) of a source with a tem- 
perature of T = 10 4 K, size R s = 2.5 x 10 13 cm (r = 0.1), 
and magnetic field B — 100 G for a — 1 (black), a = 1/3 
(red), and with l c = 10~ 2 R S and a = 1 (blue). The interac- 
tion length (thin solid), the Larmor radius (dotted), and the 
typical length-scale of energy losses from synchrotron radia- 
tion (dashed) and e + e _ pair production (dotted-dashed) in 
the case of protons are also shown. 



B = 100 G for various different diffusion regimes: The 
case of a magnetic field coherent over the entire source, 
l c = R s , is shown for Bohm diffusion, a = 1, as a black 
solid and for Kolmogorov diffusion, a = 1/3, as a red 
solid line, respectively. Additionally, i? c ff is shown for 
a = 1 and l c = 10 _2 i? s as a blue solid line and compared 
to other typical length-scale as the interaction length, 
the Larmor radius (dotted), and the synchrotron losses 
(dashed) for protons. The behavior of i? c fr is the expected 
one: Protons with E > El ph 8 x 10 17 eV propagate 
approximately on straight-lines, and thus R c s — R s . At 
lower energies, protons start to diffuse and R c g increases. 
The steepness of this increase depends on the energy de- 
pendence of the diffusion coefficient D(E), i.e. on a and 
the value of the coherence length l c . For Rl(E) < l c , 
the effective size is proportional to E~ a . If the coherence 
length turns out to be smaller than i? s , then i? c ff is pro- 
portional to E~~ 2 in the energy range from El down to 
the energy where Rl(E) = l c . In this case the effect of 
diffusion is strongest, see Fig. [5] 

The increase of R a R implies a change in the final proton 
flux, only if El is larger than the threshold energy E t h 
for photo-pion production. In Fig. we show the width 
A diff = log(£ , L )-log(-B t h) assuming E th = m p m 7T /(2e J ), 
with e 7 = 2.7T. Inside this window, even a transpar- 
ent source can become opaque, as it is illustrated by the 
increase of the effective interaction depth T e s shown in 
Fig. [5] in the case of diffusion. 
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FIG. 6: Effective interaction depth r c ff for the same source 
parameters as in Fig. [5] 



The impact of a change of T e s on the final nucleon 
fluxes is shown in Fig. [7l where unnormalizcd fluxes of 
the initial and final protons as well as of neutrons escap- 
ing from the source before decaying are shown for the 
same cases as above. The main effect of diffusion on the 
final proton spectrum is a reduction of the flux in the 
energy window between E t ^ and El- This reduction is 
proportional to the increase of the effective interaction 
depth as long as t c h < 1. Thus the dependence of the 
proton flux can be deduced directly from the behavior of 
T e ff, which in turn depends on i? c ft- 

In the case of neutrons, changes in the flux are also 
tightly correlated to variations of i? c ff. The increase of 
r e ff for protons implies more scatterings and therefore a 
larger number of produced neutrons. Since the interac- 
tion depth of neutrons is not affected by diffusion, all 
neutrons escape from the source before decaying, if the 
source is transparent. This increased production of neu- 
trons due to the diffusion of protons stops at El. At 
intermediate energies there is a flat plateau of the neu- 
tron flux. Its slope is approximately the one of the initial 
proton flux, only slightly modified by the weak energy 
dependence of the cross section. 
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FIG. 7: Unnormalized fluxes of initial protons (dashed), neu- 
trons escaping the source before decaying (dotted), and the 
cosmic ray flux (sum of escaped protons plus protons from 
neutron decays, solid line) for the same source parameters as 
in Fig. 



this limiting energy 3 follows as 



-^syn 



A ( 

2a V B 



1 



(19) 



Above this energy, the main component of the nucleon 
flux are neutrons produced in scatterings. The energy 
losses of protons lead to a steepening of the proton flux 
and therefore also of the neutron flux. This steepening 
can be estimated for a transparent source by noting that 
the energy of protons decreases with distance I in the 
classical limit as 



E 



y = 



f 



i 



E, 1 + b Ei, 



(20) 



where b = W c \/E 2 = (2/3)a(B / B cr ) 2 . Hence, the en- 
ergy of a proton that interacts and produces a neutron 
is reduced on average by the factor (y) = J^dl y oc 1 / E- t 
(neglecting logarithmic corrections). Therefore the ex- 
ponent of the neutron spectrum is increased by one, 
(±n = Oi g + 1, compared to the generation spectrum of 
protons. 

Perhaps surprisingly, the final proton flux of a magne- 
tized source is not just a broken power-law, with a = a„ 



B. High-energy protons and energy losses 



At high energies, when the synchrotron loss length l syn 
becomes smaller than the size of the source R s , cf. Fig.0 
the proton flux is strongly suppressed. From Eq. (|10[) . 



The limiting energy E syn should be distinguished from the max- 
imal acceleration energy -E max , defined by setting synchrotron 
energy losses equal to the energy gain in an electromagnetic field, 
-Emax = [(3m 4 )/(2a 3 S)] 1 / 2 . The latter is evaluated for the mag- 
netic field in the acceleration region. 
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below and a — a g + 1 above E syn . Instead, there is 
additionally a step in the flux at E 8yn , where, e.g., the 
flux shown in Fig. [7] decreases suddenly by two orders 
of magnitude. The reason for this step is that the fi- 
nal proton flux at E < E syn consists both of escaping 
"direct" protons and of the ones produced in decays of 
escaping neutrons, whereas at E > E syn only the sec- 
ond contribution is significant. The difference in the 
relative size of the final cosmic ray flux at E > E syn 
and E < E syn corresponds to the ratio R of the neutron 
and proton flux in the plateau region, here at 10 18 eV- 
10 20 eV. It can be estimated for transparent sources as 
i?o ~ T x BR(p -> n) x (1 - (y(a g ))), where 1 - (y(a g )} 
is the spectrally averaged energy transfer from protons 
to neutrons. If the interaction depth is large, protons 
produce neutrons before losing too much energy due to 
synchrotron losses. Therefore, the step at E syn , where 
the cosmic ray flux leaving the source changes from neu- 
tron to proton dominated, increases for smaller r . This 
step in the emitted cosmic ray flux at -E S yn acts in the 
case of rather transparent sources, r « 1, effectively 
as the maximal energy of the source. The distribution of 
maximal magnetic field strengths in different sources and 
the resulting distribution of maximal energies dn/dE max 
leads thus to a steepening of the observed diffuse cosmic 
ray flux and may thereby reconcile the observed steep 
spectrum a g ~ 2.6 with the expectation a g « 2.0-2.2 
from Fermi shock acceleration 1301] . 
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FIG. 8: Unnormalized fluxes of initial protons (dashed), neu- 
trons escaping the source before decaying (dotted), and the 
cosmic ray flux (sum of escaped protons plus protons from 
neutron decays, solid line) for a source with a temperature of 
T = 10 4 K, size R s = 2.5 x 10 13 cm (r = 0.1), and magnetic 
field B = 10 G (black), 100 G (red), and 10 4 G (blue). 



C. Dependence on the magnetic field strength 

Let us discuss now the impact of the magnetic field 
strength on the nuclcon flux. Increasing B has two main 
consequences: On the one hand, the energy range Adiff 
with diffusion increases, because El oc B. On the other 
hand, synchrotron radiation losses become important for 
lower energies, as E syn oc 1/B 2 . Thus an increase of 
the magnetic field strength widens the energy range Adiff 
with diffusion until El = E syn . Increasing B even fur- 
ther, Adiff narrows again and synchrotron losses influence 
finally the whole spectrum. This dependence is clearly il- 
lustrated by Fig. [51 where the nucleon fluxes for a source 
with a temperature of T = 10 4 K, size R s = 2.5 x 10 13 cm 
(r = 0.1), and three different magnetic fields B = 10 G 
(black), 100 G (red), and 10 4 G (blue) are shown. In the 
last case, the cosmic ray flux is strongly suppressed and 
the main part of the source luminosity is damped into 
electromagnetic cascades. 
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D. Dependence on the interaction depth 

Let us next consider the consequences of the magnetic 
field strength on the nucleon flux from a thick source, 
To > 1. As discussed in Ref. Q, the characteristic prop- 
erty of a thick source is the strong suppression of the 
nucleon flux at energies above the photopion production 



FIG. 9: Unnormalized fluxes of initial protons (dashed), neu- 
trons escaping the source before decaying (dotted), and the 
cosmic ray flux (sum of escaped protons plus protons from 
neutron decays, solid line) for a source with a temperature of 
T = 10 4 K, size R s = 2.5 x 10 13 cm (r = 0.1), and magnetic 
field B = 10 2 G (black), R s = 2.5 x 10 16 cm (r = 10 2 ) and 
B = 10~ 5 G (red), and R s = 2.5 x 10 16 cm (r = 10 2 ) and 
B — 10 2 G (blue). 
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threshold due to multiple scattering. For the unmagnc- 
tized, thick source in Fig. [51 the final proton flux is cut 
off at £th ~ 10 16 eV followed by a small bump, while at 
lower energies the proton flux is hardly affected as long 
as R s is smaller than the typical length scale for e + e~ 
pair production, see Fig. [5l The neutron flux is narrowly 
peaked around E t h, where the neutron interaction length 
decreases and the escape probability increases. 

The consequences of magnetic fields in a thick source 
are twofold. At high energies, there is a competition 
between multiple scattering and synchrotron radiation 
as the main energy loss process. Analogous to Eq. (fH)]) . 
we define as critical energy 



3 / B CI 
2^ \~B 



1 

hnt 



(21) 



While both processes lead to a strong suppression of the 
nucleon flux at high energies, only the latter transfers the 
energy to neutrinos. In contrast, magnetic fields influence 
not only the neutrino but also the cosmic ray flux at 
low energies. Since diffusion increases the effective size 
i? e fi of the source below El, synchrotron and e + e~ pair 
production losses can suppress the proton flux even below 
the pion-production threshold. As result, the final cosmic 
ray flux may consist only of neutrons which escape near 
E t h- As an example for such a situation, we show the 
proton and neutron flux from a source with B = 100 G 
and To = 100 in Fig. [5] that is narrowly concentrated 
around the threshold energy. 



IV. NEUTRINO YIELDS 

In order to analyze the effect of the magnetic fields 
of the source on the neutrino spectrum we study now 
the neutrino yield of a single source. The neutrino yield 
Y V (E), i.e. the ratio Y„(E) = ^(E)/( Pint (f> p (E)) of the 
emitted neutrino flux (f> u and the product of the interac- 
tion probability p; nt = I — exp(— To) and the injected pro- 
ton flux <j) p , represents the number of neutrinos produced 
per injected interacting proton with the same energy 4 . 



A. General characteristics 

The dependence of the neutrino yields on the source 
parameters can be deduced from the corresponding be- 
havior of secondary mesons in the source. Therefore 
we complement the typical length-scales shown in Fig. [5] 
for protons with those for mesons in Fig. [TU1 We have 
added the typical length of inverse Compton scattering 
for muons as well as the decay lengths l^cc = lh/2 of 
mesons and muons. 



4 Note that this definition differs from the standard one, Y U (E) 
<j) v {E) / (TQ(pp(E)), that was also used in Ref. [jj. 
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FIG. 10: Decay lengths (thick solid) and the typical length- 
scale of synchrotron losses (dashed) for charged pions (black) , 
charged kaons (blue) and muons (red) for a source with tem- 
perature T = 10 4 K, size R s = 2.5 x 10 13 cm (r = 0.1), and 
magnetic field B = 10 2 G. The interaction length for pions 
(thin solid) as well as the typical length of inverse Compton 
scattering for muons (dotted-dashed) have been included. 



The most important difference with respect to a source 
with negligible magnetic fields is the appearance of a 
length scale, Z syn cx 1/E, which decreases with energy (as 
long as x <C 1). Hence, the synchrotron length becomes 
for a magnetized source at high enough energies always 
the smallest length scale, and in turn meson decays and 
the neutrino flux will be suppressed. 

We modify Eq. (fl9)) to take into account possible me- 
son and muon decays, defining 



Ei. 




for fi(B, R s ) <l 



for fi(B, R s 



where the functions fi(B,R s ) are defined as 



fi(B,R s ) 



2a 

T 



B 

Bin 



Rim, 



I, 



1/2 



>1, 

(22) 
(23) 



The first expression represents the condition l S y n (E syn ) = 
R s . At energies above E syn , charged mesons lose energy 
by synchrotron radiation. Mesons with energies between 
Rslijilm and -E S yn escape from the source before decay- 
ing. At energies smaller than Rg^^/m mesons decay 
inside the source. The second expression corresponds to 
l S yn{E S yn) = 'dcc(£<syn)- In this situation mesons with 
energies below E syn decay always within the source. 

For given source parameters, the limiting energy E syn 
above which synchrotron losses become crucial increases 
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FIG. 11: Neutrino yield Y v from the decay of charged pions 
(solid red), charged kaons (solid blue), if£ (dashed blue), K\ 
(dotted blue), and neutrons (solid magenta). In the case pions 
we also show the fraction of neutrinos produced by the muons 
generated in the pion decay (dotted red). The same source 
parameters as in Fig. 1101 have been assumed. 

with mass and decreases with the decay length of a 
particle. As a consequence, muons and pions are af- 
fected already at much lower energies than kaons and 
protons [IlGj- 

In Fig. [Til we show how different reactions contribut- 
ing to the neutrino yield are affected by magnetic fields. 
The increase of the number of ^7 interactions at low ener- 
gies by diffusion has two main consequences for the neu- 
trino production: First, a larger number of neutrons is 
produced, which will freely escape in the case of a trans- 
parent source. Second, the number of generated pions 
increases. The former leads to a bump of v e around the 
energy f» 10 _3 i?th- The latter produces another bump of 
neutrinos produced by pion decay roughly at E t h- These 
two characteristics can be clearly observed in Fig. [TT1 

The high energy end of the spectrum is mainly influ- 
enced by synchrotron losses. The relative importance 
of different channels to the total neutrino yield is de- 
termined by the relative size of the synchrotron energy, 
Eg yn < Es yn < E^ D , cf. Fig. [Til] Therefore, between 
the threshold energy for pion production and E£ yn pion 
decays are the most important contribution to the neu- 
trino yield. At slightly higher energies, between E£ yn and 
E^ n , neutrinos from the decay of charged kaons domi- 
nate. For the source considered in Fig. [ITjJ all neutrinos 
are produced in the decay of A£'s at the highest energies, 
as they are not directly affected by synchrotron losses. As 
discussed for neutrons, synchrotron losses of protons lead 
however to a steepening of the A'£ spectrum and in turn 
also to a steeper neutrino spectrum from these decays. 



We note that, in contrast to K\, the K s contribution is 
suppressed roughly at E£ yn . Reason for this is that the 
main decay channel of Kg's giving rise to neutrinos is 
ir + ir~ . 

Beside the neutron and pion neutrino bump at low 
energies and the suppression at high energies due to syn- 
chrotron losses, there is a third characteristic feature 
at intermediate energies: Muon damping in the energy 
range between Eg yn and E£ yn [lj| . While the synchrotron 
length l syn of pions and muons are comparable, their de- 
cay lengths /doc differ by more than two orders of magni- 
tude. As a consequence, muons with an energy between 
£^, n and E£ yn created from pion decay can not further 
decay but lose energy by synchrotron radiation. Hence, 
most neutrinos in that energy range are produced di- 
rectly in pion decays. 

In order to estimate the energy range where the ef- 
fect of muon damping is potentially strong, we define the 
muon damping width A M analogous to Adiff as 

A M EE log(A- n ) - log(^yn) • (24) 

The muon damping width depends on the source param- 
eters R s and B only indirectly via the choice between 
(p?2a ) and (|22b ). Let us first consider a source such that 
/ M (B,i? s ) > 1, i.c (B/G)(R s /cm) > 3 x 10 15 . In this 
case, all muons produced in pion decay with energies be- 
tween Eg and E£ are generated inside the source and 
therefore are affected by the synchrotron losses. Hence, 
muons are completely damped within this energy range. 
The muon damping width is maximal and equals A M ss 
0.5 x log^/m^) 5 ^/^)] » 1.3. If U(B,R S ) < 1, 
some of the muons produced can escape from the source 
before decaying. Let us now consider a source such that 
U(B,R a ) < 1, i.c. (B/G)(R s /cm) « 5x 10 14 . In 
this case, pions with energies between Rgl^fi/m^ and 
E£ yn decay outside the source, and then there is no 
muon damping. The numerical values used in Fig. 111! 
(B/G)(R s /cm) = 2.5 x 10 15 , correspond to partial muon 
damping, as can be observed in the figure. 

A similar damping effect occurs for muons from 
charged kaon decays. Now the numerical value of the 
width Eq. ([24]) is larger. However, the effect is not as 
strong as in the case of pions, because the energy range 
where kaons dominate is narrower and also their number 
of decay channels is larger. 

B. Dependence on the properties of the source 

Wc now analyze how the general characteristics previ- 
ously discussed depend on the different properties of the 
source. 

1. Diffusion regime 

In Fig. 1121 wc present the neutrino yields for sources 
with different diffusion conditions, choosing those already 



11 




FIG. 12: Neutrino yield Y v produced in sources with the same 
parameters considered in Fig. [5] The different reactions giving 
rise to Y v in the case of the source with Bohm diffusion (a = 1) 
and l c = R 3 are also shown. 

used in Scc. lIII Al As previously discussed, diffusion gen- 
erates two bumps, one from neutrinos created by neutron 
decay at very low energies, and another one of neutrinos 
produced in the decay of charged pions at E m Eth . The 
intensity of these peaks is directly related to the charac- 
teristics of diffusion, in particular to r c ff. A comparison 
of the yields with Figs. [5] and [6] confirms the expected cor- 
relation between the height of the neutrino bumps at low 
energies and the effective interaction depth at the cor- 
rect energy range. Also, the different diffusion regimes 
do lead to the same neutrino yield at high energies. 



FIG. 13: Neutrino yield Y v produced in sources with the same 
parameters assumed in Fig. [8] The different contributions to 
Y v in the case of the source with B = 10 G and B — 10 4 G are 
also shown. The fraction of neutrinos produced by the muons 
generated in the pion decay (dotted red) is also indicated. 



suppressed in the whole energy range, cf. Fig. 1131 

The value of the magnetic field strength also affects 
muon damping, enhancing the damping for a stronger 
field. For the case of a source with 10 G shown in Fig. Q21 
the product (B/G)(R s /cm) = 2.5 x 10 14 < 5 x 10 14 . 
Therefore most pions decay outside the source and thus 
muons arc not affected by the magnetic field. If the mag- 
netic field is increased to 10 4 G, then (B/G)(R s /cm) = 
2.5 x 10 17 ^> 3x 10 15 and complete muon damping occurs. 



2. Magnetic field strength 

In Fig. 1131 we show the dependence of the neutrino 
yield on the magnetic field strength. As discussed in 
Sec. lIII Cl two main effects are expected. At low energies, 
the larger the magnetic field strength the more intense is 
diffusion and hence the interaction depth. This translates 
into an increase of the neutrino yield at low energies, and 
in particular into a higher bump of the contribution to the 
neutrino yield coming from pion and neutron decay. This 
can be easily confirmed by comparing the two sources 
with 10 G and 100 G in Fig. [H 

The second consequence of a larger magnetic field 
strength is a more severe suppression of the neutrino yield 
at high energies due to synchrotron radiation. Since the 
energy at which synchrotron losses become important is 
proportional to 1/-B 2 , the suppression of the neutrino 
yield starts earlier in the source with 10 2 G than in the 
one with 10 G, while for B = 10 4 G the neutrino yield is 



3. Interaction depth 

We analyze now the influence of the interaction depth 
t on a magnetized neutrino source. Figure [T4l shows the 
different contributions to the neutrino yield for a source 
with the same parameters as the one in Fig. QTJ but for 
t = 100 instead of t = 0.1. 

At low energies, the narrowly peaked spectrum of neu- 
trons leaving a thick source is responsible for the absence 
of the high energy tail of the neutrinos produced in neu- 
tron decay. The number of neutrinos from pion decay 
in this range is not influenced strongly, because multiple 
scattering occurs in both cases either because of large T e ff 
or large To. 

At high energies there are some differences. In the 
case of thin sources the most important contribution to 
the neutrino yield comes from the decay of neutral kaons, 
see Fig. [TTJ since they are not directly affected by syn- 
chrotron losses. In the case of sufficiently thick sources, 
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FIG. 14: Neutrino yield Y v from the decay of charged pions 
(solid red), charged kaons (solid blue), if£ (dashed blue), K\ 
(dotted blue), neutrons (solid magenta), and charm mesons 
(solid green) for a source with temperature T = 10 4 K, size 
R s = 2.5 x 10 16 cm (r = 10 2 ), and magnetic field B = 10 2 G. 
The fraction of neutrinos produced by muons from the pion 
decay (dotted red) is also indicated. 



also most neutral kaons interact with the photons and 
thus the neutrino production from these decays is sup- 
pressed, see Fig. [9l In Ref. Q, it was shown that neu- 
trinos from decays of charmed mesons became the domi- 
nant contribution at the highest energies for a source with 
To 3> 1. In principle, the same holds in the case of sources 
with magnetic fields. However, synchrotron losses above 
E > -Esyn may significantly suppress the proton flux and 
thus also the production of charmed mesons produced. 
For the source considered in Fig. [51 the contribution from 
charm decays is negligible. 



In the case of charged mesons and muons it is possible 
to redefine the energy at which synchrotron radiation be- 
comes the dominant energy loss process substituting R s 
by iint into Eq. (f2"2")) . If Z;„t) < 1 then pions with en- 

ergies between knt^n/lo and E^ yn scatter before decay- 
ing. If / T (B,iint) ^ 1 then pions decay inside the source 
whereas muons with energies higher than E£ yn will suffer 
energy losses by synchrotron radiation before decaying. 
Therefore the condition for having muon damping in a 
thick source can be written as f- K {B 1 li nt ) > 1 or cquiva- 
lently (B/G) (l- mt /cm) >5x 10 14 , which clearly satisfied 
by the source considered in Fig. [9] 



V. FLAVOR DEPENDENCE OF NEUTRINO 
YIELDS 

Both neutrino telescopes and extensive air shower ex- 
periments have some flavor discrimination possibilities. 
The case of neutrino telescopes was discussed for the ex- 
ample of ICECUBE in detail in Ref. [3l[ , for air-shower 
experiments see Ref. (32|. The main observable for the 
neutrino flavor composition is the ratio of track to shower 
events in a neutrino telescope, = (j)^/ (</> e + <fi T ), while 
only in a very small energy range all flavors can be distin- 
guished. Additionally, extensive air shower experiments 
are sensitive to the fraction of tau events in all horizontal 
neutrino events, R T = (f> T /((f> e + in a small energy 
window around 10 18 eV. 

In spite of these flavor discrimination possibilities, the 
potential of high energy neutrino observations for mixing 
parameter studies has been realized only recently (33j , for 
some follow-up studies see [IH ■ One of the main reasons 
for this missing interest has been the prejudice that the 
maximal mu-tau mixing together with the expected fla- 
vor ratio <j){v e ) '■ 4>{ v ^) '■ 4>( v t) = 1:2:0 from pion 
decay prevents oscillation studies with high energy neu- 
trinos. Key observation of Refs. [33| was that there exist 
however several examples of neutrino sources where at 
least in some energy range significant deviations from this 
canonical flavor ratio from pion decay can be expected. 
Our results indicate that deviations from the canonical 
flavor ratio <f>(v e ) '■ 4>( v n) '■ ^(^v) = 1:2:0 are much 
more common than previously thought. In the following, 
we will show that magnetized sources are characterized 
by a strongly energy-dependent flavor ratio and thus the 
flavor ration encodes non-trivial information. 



A. Energy dependence of the neutrino flavor ratio 
at the source 

We start with an analysis of the expected neutrino fla- 
vor ratio R° = Y^^/Y^^ at the source. In Fig.[15]we show 
the flavor ratio at the source Rq for four different sources 
withT= 10 4 : B = 0G, and R s = 2.5xl0 13 cm (r = 0.1) 
(black), B = G, and R s = 2.5 x 10 16 cm (r = 100) (ma- 
genta), B = 100 G, and R s = 2.5 x 10 13 cm (r = 0.1) 
(blue), B = 100 G and R s = 2.5 x 10 16 cm (r = 10 2 ) 
(red). In the case of a transparent source with negligi- 
ble magnetic field strength the flavor ratio basically does 
not depend on energy. However, the presence of mag- 
netic fields as well as large interaction depths can result 
in complicated structures of the flavor ratio. The differ- 
ent bumps arising at different energies are tightly related 
to the different processes taking place in the formation 
of the neutrino flux. 

Let us analyze these structures for each case. The 
flavor ratio Rq of all sources is zero below the pion- 
production threshold as most neutrinos are v e produced 
in neutron decay. Right after the energy threshold for 
pion production, Rq becomes two, typical for tt decay. At 
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FIG. 15: Flavor ratio 7?° of the neutrinos produced in a source 
at T = 10 4 , B = G, and R s = 2.5 x 10 13 cm (r = 0.1) 
(black), B = G, and R s = 2.5 x 10 16 cm (r = 100) (ma- 
genta), B = 100 G, and R s = 2.5 x 10 13 cm (r = 0.1) (blue), 
and R s = 2.5 x 10 16 cm (r = 100) (red) before they oscillate. 
The main reactions giving rise to the neutrino signal are also 
shown. The neutrinos produced in the muon damping are 
indicated as direct n decay. 



higher energies the behavior of the flavor ratio strongly 
depends on the characteristics of the source. Let us first 
summarize the two cases with negligible magnetic fields. 
If the source is transparent, then Rq remains constant, 
Ro ~ 2. Thick sources show, though, a bump exactly at 
the energy at which there is a crossover between pion 
and charged kaon neutrino dominance, namely, where 
li nt (E) = l^ cc (E). At higher energies, the flavor ratio 
tends to one, when all neutrinos come from the decay of 
charmed mesons, see Ref. @. 

In the case that sources have non-negligible magnetic 
fields the energy dependence of the flavor ratio becomes 
more complex. The main feature is the presence of a 
strong peak arising between E^ yn and E£ yn , in this case 
around E f=a 10 17 eV. This clearly reflects the energy 
window where muon damping occurs. In this case the 
only neutrinos produced are the generated directly 
from pion decay. The accompanying muons lose energy 
instead of decaying, therefore no electron neutrinos are 
produced. As a consequence Rq increases, the final value 
depending on the strength of the muon damping. By 
comparing Figs. [11] and [TJJ one expects a stronger ef- 
fect for a thick source, as can be observed in the fig- 
ure. At higher energies there is an energy range where 
neutrinos are produced in the decay of charged kaons, 
analogously to the case of thick source with negligible 
magnetic fields. In contrast to this case, though, the en- 
ergy range of charged kaon dominance happens at lower 



energies. The reason is that the crossover between pions 
and charged kaon neutrinos takes place at E£ yn , which, 
for the sources considered in Fig. [T51 is smaller than the 
energy where l[ nt {E) = l^ cc (E). The different height of 
the bump observed in the transparent and thick sources 
is related to the muon damping happening also for the 
charged kaons in the latter case. Finally at the highest 
energies one expects a tendency toward one, either be- 
cause these neutrinos are produced via K® — > ir + e + v e , 
as in the figure, or because they come from the decay of 
charmed mesons. 

We conclude that there is a direct relationship between 
the energy dependence of the flavor ratio Rq and the dif- 
ferent processes contributing to the neutrino yield, that 
may be used to deduce the intrinsic characteristics of the 
source. 



B. Effect of neutrino oscillations 

The neutrino spectra at the source discussed in the pre- 
ceding subsection are modulated by oscillations. There- 
fore the expected flavor ratios Ri at the Earth are differ- 
ent from the original ratios at the source. 

The neutrino fluxes arriving at the detector, <f>® , can 
be written in terms of the initial fluxes <p a and the con- 
version probabilities P a p, 







1 ae 



Ear 



(25) 



Since the interference terms sensitive to the mass split- 
tings Am 2 's do not contribute, the conversion probabili- 
ties are simply 



p a0 = s a p-2j2 mh u 0kU aj u* k ) 
j>k 



(26) 



where U is the neutrino mixing matrix and Greek (Latin) 
letters are used as flavor (mass) indices. 

In Fig. [TBI we show the expected flavor ratio at the 
detector for a source with T = 10 4 , B = 10 2 G, and R s = 
2.5 x 10 16 cm (t = 10 2 ), as a solid, red line for the best- 
fit point of the neutrino mixing parameters, sin 2 $12 = 
0.3, sin 2 $23 = 0.5 and $13 = [35|. Various bands 
show the range of i? M allowed if the mixing parameters 
are varied within their 95% C.L., while the unconstrained 
CP phase i5cp is varied in the whole possible range, <5cp € 
[0 : 7r] . Dominant uncertainty for the prediction of is 
the value of $ 2 3- 

The main effect of neutrino oscillations is to smear out 
the structures of the initial flavor ratio Rn. Nevertheless, 
the remaining energy dependence still allows for an iden- 
tification of the different regimes described in Sec. IV Al 
On the other hand, its dependence on the neutrino mix- 
ing parameters offers the additional possibility to obtain 
information on neutrino properties. 
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FIG. 16: Flavor ratio i? M at the Earth for a source with T = 
10 4 , B — 10 2 G, and R s = 2.5 x 10 16 cm (r = 10 2 ). The solid 
red line (gray) corresponds to the best-fit point of the neutrino 
mixing parameters, sin 2 $12 = 0.3, sin 2 $23 = 0.5 and $13 = 0. 
In the colored areas mixing angles within 95% C.L. and a 
possible non-zero value of <5cp have been considered: #13 and 
5c p in yellow (very light gray), plus #23 in blue (dark gray), 
and $12 in green (light gray). 



VI. SUMMARY AND CONCLUSIONS 

We have calculated in this work the yield of high en- 
ergy neutrinos produced in astrophysical sources for arbi- 
trary interaction depths tq and magnetic field strengths 
B. Diffusion of protons increases their path- length and 
hence also the effective interaction depth T e g. Therefore 
magnetized sources can lose their transparency, while the 
neutrino fluxes can be increased in this energy range. 
Sources with r e ff > 1 require to account for multiple scat- 
terings and for secondary meson-photon interactions. For 
the latter task, we have extended the SOPHIA model in 
a self-consistent way. 



Magnetic fields present in the surrounding of the ac- 
celeration region lead not only to diffusion of charged 
particles, but induce also synchrotron radiation as most 
important energy loss process. Synchrotron losses result 
in a strong suppression of the cosmic ray flux and hence 
also of the neutrino flux at high energies, E > E syn . 

The cosmic ray upper bound of Refs. @, [1[ relies on, 
among others, the assumption that the cosmic ray (with- 
out propagation effects) and the neutrino flux have the 
same spectral shape. The strong deviation of the es- 
caping cosmic ray and neutrino fluxes from the injected 
power-law found by us therefore undermines the basis of 
this bound even for transparent sources. 

Since the relative importance of the various channels 
contributing to the neutrino yields changes strongly as 
function of the energy, large variations exist in the neu- 
trino flavor composition emitted by a magnetized source. 
These variations are for magnetized sources even stronger 
than those found previously for sources with negligible 
magnetic fields in Ref. @- In particular, we have ex- 
amined two specific examples discussed earlier in the lit- 
erature, namely the cases where the neutrino spectrum 
is dominated by kaon decays [ll|, [l2j or influenced by 
the damping of muons produced in pion decays fl3j . In 
addition, we have pointed out the possible existence of 
an analogous muon damping in the case of charged kaon 
decay at high energies. We have analyzed the condi- 
tions required for a source to present these features in the 
neutrino yield. In particular, we have found that muon 
damping and kaon dominance may influence the neutrino 
yields in the same source, but in different energy ranges. 

The application of our simulation to concrete astro- 
physical models and the calculation of diffuse neutrino 
fluxes will be performed in paper III of this series [l6| . 
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